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By employing a twisted mass term, we compare recent results from lattice calculations of Nf = 2 dynamical 
Wilson fermions with Wilson Chiral Perturbation Theory (WChPT). The final goal is to determine some com- 
binations of Gasser-Leutwyler Low Energy Constants (LECs). A wide set of data with difi'erent lattice spacings 
(a ~ 0.2 — 0.12 fm), difi'erent gauge actions (Wilson plaquette, DBW2) and difierent quark masses (down to the 
lowest pion mass allowed by lattice artifacts and including negative quark masses) provide a strong check of the 
applicability of WChPT in this regime and the scaling behaviours in the continuum limit. 



1. Introduction 

It was recently recognized |1I2I3) that an im- 
portant obstacle for the simulation of Lattice 
QCD with light dynamical quarks comes from the 
phase structure of Lattice QCD at small quark 
masses. Lattice artifacts of O(a^) may prevent 
the pion mass from being smaller than a certain 
minimum value. Such a scenario had been pre- 
dicted - as an alternative to an Aoki phase sce- 
nario - in the context of Chiral Perturbation The- 
ory for the Wilson fermion action (WChPT) 0], 
and confirmed numerically in These possi- 

bilities are not only a feature of all Wilson-like 
fermions, but also of staggered fermions |5I6I7| . 
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The most disturbing aspect of this problem is that 
it may be easily overlooked. In fact very long 
Monte Carlo histories have been observed, which 
eventually turned out to be meta-stable points 
when jumping to a higher value of the pion mass 
In the case of pure Wilson fermions long 
metastable HMC histories (^ 1000 trajectories) 
have been recently observed (before converging 
to a single point) even for lattice spacings as low 
as a « 0.08 fm and pion masses m^r ~ 300 MeV 
|Hj. Even if dangerous, the problem has a simple 
solution (although for an additional cost): the 
comparison of simulations at positive and nega- 
tive quark masses allows to recognize where meta- 
stabilities appear. 

In this scenario the effect of the gauge action 
has proved to play an important role. In par- 
ticular simulations with the DBW2 |5] gauge ac- 
tion showed a considerable decrease of the mini- 
mal pion mass [2]. In this work we look in detail 
into this effect. 

The twisted mass fermion approach jJOl 
(tmQCD) provides, among other advantages, the 
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ideal framework for the investigation of the (zero- 
temperature) phase diagram of lattice QCD with 
Wilson fermions, see refs. |11I12I13| for reviews on 
twisted mass fermions in present and past confer- 
ences. On the analytical side WChPT (4iI4iIB] . 
which has been extended to the twisted mass case 
|ie)ll7ll8llHl^ . offers an efficient tool to inter- 
pret the lattice data. 

Recently we have collected a large statistics of 
lattice data at large and moderate lattice spacing 
(a ~ 0.2 — 0.12 fm), different gauge actions (Wil- 
son plaquette and DBW2) and different quark 
masses (down to the lowest pion mass allowed 
by lattice artifacts and including negative quark 
masses) . In this work we show explicitly the com- 
parison between these lattice data and WChPT. 
Although for reliable physical predictions defi- 
nitely smaller lattice spacings are needed, also 
the lattice spacings that we consider here are use- 
ful as a starting point for the extrapolation to a 
continuum limit. Nevertheless, the results of the 
present work have to be seen as an exploratory 
study. 

By now a good amount of data has been col- 
lected also for the tree-level Symanzik improved 
gauge action (tlSym) |21| , which are very promis- 
ing, but a detailed ChPT analysis is not yet ready. 
We will only comment briefly on that at the end. 

2. Lattice simulations on Nf = 2 QCD 

2.1. Fermionic and gauge action 

The lattice action for a doublet of degenerate 
twisted mass Wilson fermions (in the so-called 
"twisted basis") reads 

Sq = ^{{XxifJ''^ + ilbTza^i]xx) + (1) 

X 

1 

with = amo + Ar — 1/(2k), r the Wilson- 
parameter, set in our simulations to r = 1, arriQ 
the bare "untwisted" quark mass in lattice units 
(k is the conventional hopping parameter) and /i 
the twisted quark mass; we also define Ux.-fi = 

uLf^^f, and = -7m- ^ 

The first reason for including /i is that the 



fermionic determinant is free from exceptional 
configurations when /i 7^ llQj. which helps to 
perform numerical simulations with light quarks. 
The second reason is a general 0{a) improvement 
[22 and the reduced operator mixing (23i24] . 
when = /XK,crit- 

For the gauge sector we consider the one- 
parameter family of actions including planar rect- 
angular (1 X 2) Wilson loops (C/^^^): 

4 1 

+ (2) 

with the normalization condition cq — 1 — 8ci. 
We considered three cases: i.) Wilson plaquette 
gauge action, ci = 0, ii.) DBW2 gauge action |2, 
ci = —1.4088, iii.) tree-level Symanzik improved 
gauge action (tlSym) j^J, ci = —1/12. 

The reason why such variations have been ex- 
plored relies on the experience that this may 
improve the spectrum of the Wilson fermion op- 
erator. 

2.2. Analysis of the data in tmQCD 

In this section we review some formulae which 
are relevant for the analysis of lattice data pro- 
duced with tmQCD (2). In fact tmQCD offers 
new possibilities for the determinations of the ba- 
sic QCD parameters and renormalization factors. 

The twist angle w defines the chiral rotation 
relating twisted mass QCD to ordinary QCD. In 
the case of the vector and axial currents the ro- 
tation reads (considering only charged currents, 
a = 1,2): 

Vxt. = ZvV^^co^uj +eabZAAl^Bmuj, (3) 
= ZaAI^ cosuj +eabZvV^^^su^uj, (4) 

where the hatted currents on the l.h.s. denote the 
chiral currents of QCD (physical currents), while 
the currents on the r.h.s. are the corresponding 
bilinears of the quark-field in the twisted (x-) ba- 
sis. Note that the renormalization constants of 
these bilinears, Zy and Za, are involved. For a 
given choice of the lattice parameters, the twist 
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angle lu is determined by requiring parity conser- 
vation for matrix elements of the physical cur- 
rents |26I2| . Since unknown renormalization con- 
stants are involved, two conditions are required, 
our choice being: 



mass QCD we use 



X 



= 




(5) 
(6) 



The solution of Eqs. lO and © with Eqs. ^ 
and 10} gives a direct determination of the twist 
angle to and of the ratio Za/Zv from lattice data, 
see for details. In particular at full twist where 
u! — tt/2 the condition reads y~^(yl^Q P~) = 0. 

X 

The full twist situation can be also obtained by 
requiring the vanishing of the PCAC quark mass 
in the x basis. Both definitions are optimal in the 
sense of |2()l1<jl2V) . 

The knowledge of the twist angle is necessary 
for the determination of physical quantities like 
the quark mass and the pion decay constant. The 
physical PCAC quark mass m^'~^^'~^ can be ob- 
tained from the Ward identity for the physical 
axial- vector current. An interesting possibility is 
to use Eqs. and parity restoration together 
with the conserved vector current of the x-fields 
1. This gives 



for which Zv 



^PCAC 



1 (V*F+P~ 
''2smuj {P+Py) 



(7) 



Analogously, for the physical pion decay constant 
/tt we use 



a^ = (am,)-l(0|i+(0)|7^^ 



(0|K)+(0)k+) ^ 
[aniT:) sin a; 

(8) 



Notice that with this definition the lattice de- 
termination of has automatically the correct 
normalization 28 29 . 

Finally, the renormalization constant of the 
vector current Zy can be determined on the basis 
of the non-renormalization property of the con- 
served current Vx^ [301 ■ In the case of twisted 



Zv = 



(0|V^+|^+) 



(9) 



We also observe that in a mass independent 
renormalization scheme, the renormalization fac- 
tors need to be extrapolated to the chiral limit 
(m^*^'^'-^ — > 0) . For this choice of renormalization 
the parity conserving conditions © , 10 , hold up 
to lattice artifacts. 

In this work we consider the lattice data with 
parameters summarized in Tabled The scaling 
behaviour of the data has already been shown 
in |31| . Here we concentrate on the comparison 
with ChPT, which is discussed in the next sec- 
tion. 



Table 1 

Simulation points. 



Action 


/3 


a [fm] 




L '"^ xT 


DBW2 


0.67 


0.19 


0.01 


12^ X 24 


DBW2 


0.74 


0.12 


0.0075 


16^ X 32 


plaq 


5.1 


0.20 


0.013 


12^ X 24 


plaq 


5.2 


0.16 


0.01 


12^ X 24 


plaq 


5.3 


0.14 


0.008 


16^ X 32 



3. WChPT versus LQCD 

The extension of WChPT to the case of adding 
a twisted mass term was considered in refs. [321161 
I17ll9llkl!j0| . We denote by L, the usual Gasser- 
Leutwyler coefficients p^. while W, W and W 
are some combinations of the new LECs as- 
sociated to 0{ma) and 0{a?) lattice artifacts 

The new LECs depend on the lattice action and 
also - in general - on the definition of the mass 
parameter, and of the currents. 

If 

Kcrit is chosen - for instance - from the van- 
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ishing of cos(w) and if we define (/i/j — ii/Zp) 



V XH+' « 

then we have for the pion mass and the PCAC 
quark mass: 



1 



In 



(47r^^o) 



+ 



--2{(4L6 + 2L8-2L4-i5)x' 



+2(2VF- W^)pXCOS {u 
+AW'p^ cos (cj)^}, 



,PCAC 



Zp 

+2W'p^ cos {ujf)] 



[X+ -^{WxP cos (w) + 
^0 



(10) 
(11) 



Similar formulae are available for /^r and (/tt 
These formulae are valid in the regime where 
m^^^^ / Kqcd ^ o-Aqcd, where most of our sim- 
ulated points are located. In the case of a first or- 
der phase transition scenario (as it appears to be 
the case) , the same formulae hold also for smaller 
masses (although some terms can be dropped). In 
the regime of very small masses also NLO calcu- 
lations have been done, leading to the addition of 
the 0{a^) terms, [34135) . however we will neglect 
these corrections here. 

We first show a qualitative comparison of our 
data with the formulae above. Figure (left) 
displays the pion mass in the case of DBW2 
gauge action and a relatively large lattice spacing 
(a ~ 0.19 fm). The data for can be reason- 
ably fitted by straight lines (therefore neglecting 
NLO ChPT terms and the small twisted mass /i) . 
What cannot be neglected is the presence of a 
minimal pion mass of roughly ~ 300 MeV. No 
meta-stabilities have been detected. If we com- 
pare with Figure n (right) - whose data are pro- 
duced with plaquette gauge action at an even 
smaller lattice spacing (a ^ 0.16 fm) - the strik- 
ing difference is the presence of a much higher 
minimal pion mass. Here many meta-stable 



points have been detected, which did not tunnel 
into stable ones. Notice that the comparison is 
done at the same value of fi, and a/i is larger for 
the DBW2 case (the effect of ^ is however small, 
at the present value). 

This picture is confirmed by the PCAC quark 
mass. Eq. (111(1 shows that lattice artifacts break 
the linear relation between m^'-^^'~^ and x (ihsy 
are simply two definitions of the quark mass). In 
particular the W term induces a difference in the 
slope at positive and negative quark masses. The 
term W', instead, is responsible for a jump at 
the origin, which prevents \m^'^^'^\ from being 
smaller than some 0{a^) lattice artifact. Both 
these expected features are visible, but small, in 
Figure 13 (left): two different slopes and a small 
jump. Also FigureEl (right) can be understood in 
terms of such effects, which however completely 
upset the continuum picture. The scaling be- 
haviour of these effects is confirmed by the other 
(3 values in Tabled 

The observation that the optimal critical mass 
can be determined by the vanishing of to^'^"^'^ 
suggested to express the pion mass (and in 
general all pionic quantities) as function of the 
PCAC quark mass instead of 1/k. The result for 

TOtt is (now X '■— '^'^ \J {cos{uj)m^^ Y + ^j?, 

other definitions are as above): 

1 



X 



32tt^F^ 



X In 



X 



(4^Fo 



+ -2 {(4^6 + 2Ls - 2U - L^)X^ 

+2{W -W)px cos (uj)}, (12) 

In this reparametrization the constant W' dis- 
appears, - and so does the 0(0^) term - and 
the pion mass can apparently go to zero when 
^PCAC _^ Q However, one should keep in mind 
that not all values of m^'~'^'-^ are accessible with 
stable simulation points. This parametrization 
allows to include in the ChPT fit also meta-stable 
points, where both m.^ and nig'^^'^' are lower 
than it would be possible in a stable minimum 
of the effective potential. Since this is an inter- 
esting check, we exploit this possibility and we 
include also meta-stable points (from |3|) in the 
fit. 
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Figure 1. Pion mass for DBW2 (left) and plaquette (right) gauge action at /3 = 0.67 (DBW2), and 
(3 = 5.2 (plaquette). For both V= 12^ x 24 and afi = 0.01. Fit of is done with a linear fit. 



Combined fits of m,r, /n and g,^ provide strong 
constraints. Results for are shown in Fig- 
ure Ol They give an estimate of Fq ~ 85 MeV 
and values for the LECs in agreement with previ- 
ous estimates [36137) . Results for the W's appear 
compatible with zero. Details will be presented 
elsewhere [55] . 

4. Conclusions 

We have shown a comparison of unquenched 
lattice data with WChPT. We found that 
WChPT seems to describe the data rather well 
and results for low-energy constants are consis- 
tent with previous determinations. This lets us 
be confident that a precise determination will be 
possible with reasonable computation cost. 

WChPT lets us investigate the phase structure 
of lattice QCD. There is a striking difference on 
this aspect between DBW2 and plaquette gauge 
action. It is interesting to vary the coupling ci 
in Eq. (0) and interpolate between DBW2 and 
plaquette actions. It appears that even a small 
value of ci can already have a large impact on 
the phase structure. This is illustrated in Fig- 
ure 01 where we show the average plaquette value 
as a function of the hopping parameter k for 



three different actions, i.e. different values of ci, 
namely ci = (Wilson), ci = —1/12 (tlSym) and 
ci = -1.4088 (DBW2). As one moves k from the 




0.165 0.17 0.175 0.18 



Figure 4. Hysteresis of the average plaquette 
value as n is moved across the critical point, 
for Wilson, tlSym and DBW2 gauge action at 
a - 0.17 fm. 



negative or positive side across the critical point. 
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Figure 2. PCAC quark mass for DBW2 (left) and plaquette (right) gauge action at /3 = 0.67 (DBW2), 
and /3 = 5.2 (plaquette). For both V= 12^ x 24 and a/i = 0.01. Fits are only qualitative, on the basis of 
Eq. 



where the PCAC quark mass vanishes, a hystere- 
sis in the average plaquette value develops whose 
size and width are indicators of the strength of 
the phase transition. We observe that both the 
width and the size of the gap in the plaquette 
value decreases considerably as we switch on ci 
to ci = — 1/12 (tlSym action). Decreasing ci fur- 
ther down to ci = -1.4088 (DBW2 action) still 
seems to reduce the size of the gap, but the effect 
is surprisingly small despite the large change in 
ci . Note that the results in Figure 01 are for a 
lattice spacing a ~ 0.17 fm that is roughly con- 
sistent for all three actions. One should remark 
that the results for the Wilson plaquette gauge 
action are at non-zero twisted mass fi = 0.013, as 
opposed to the tlSym and the DBW2 data in the 
same plot. Since the strength of the phase tran- 
sition is expected to be reduced as one switches 
to a non-zero twisted mass, a true comparison at 
/i = would disfavor the Wilson plaquette gauge 
action even more. 

Comparison of tlSym data with ChPT will be 
presented elsewhere IBP]. 
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